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Abstract — Cellular networks are in a major transition from 
a carefully planned set of large tower-mounted base-stations 
(BSs) to an irregular deployment of heterogeneous infrastruc- 
ture elements that often additionally includes micro, pico, and 
femtocells, as well as distributed antennas. In this paper, we 
develop a tractable, flexible, and accurate model for a downlink 
heterogeneous cellular network (HCN) consisting of K tiers of 
randomly located BSs, where each tier may differ in terms of 
average transmit power, supported data rate and BS density. 
Assuming a mobile user connects to the strongest candidate 
BS, the resulting Signal-to-Interference-plus-Noise-Ratio (SINR) 
is greater than 1 when in coverage, Rayleigh fading, we derive an 
expression for the probability of coverage (equivalently outage) 
over the entire network under both open and closed access, which 
assumes a strikingly simple closed-form in the high SINR regime 
and is accurate down to — 4 dB even under weaker assumptions. 
For external validation, we compare against an actual LTE 
network (for tier 1) with the other K — 1 tiers being modeled 
as independent Poisson Point Processes. In this case as well, our 
model is accurate to within 1-2 dB. We also derive the average 
rate achieved by a randomly located mobile and the average load 
on each tier of BSs. One interesting observation for interference- 
limited open access networks is that at a given SINR, adding more 
tiers and/or BSs neither increases nor decreases the probability of 
coverage or outage when all the tiers have the same target-SINR. 

Index Terms — Femtocells, heterogeneous cellular networks, 
stochastic geometry, point process theory, coverage probability. 

I. Introduction 

Mathematical analysis for conventional (1-tier) cellular net- 
works is known to be hard, and so highly simplified system 
models or complex system level simulations are generally 
used for analysis and design, respectively. To make matters 
worse, cellular networks are becoming increasingly complex 
due to the deployment of multiple classes of BSs that have 
distinctly different traits (3). For example, a typical 3G 
or 4G cellular network already has traditional BSs that are 
long-range and guarantee near-universal coverage; operator- 
managed picocells Q, |5) and distributed antennas (6|-||9} 
that have a more compact form factor, a smaller coverage area, 
and are used to increase capacity while eliminating coverage 
deadzones; and femtocells, which have emerged more recently 
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and are distinguished by their end-user installation in arbitrary 
locations, very short range, and possibility of having a closed- 
subscriber group (T0|-| 12 1. This evolution toward heterogene- 



ity will continue to accelerate due to crushing demands for 
mobile data traffic caused by the proliferation of data-hungry 
devices and applications. 




Fig. 1. Illustration of a three-tier heterogenous network utilizing a mix of 
macro, pico and femtocell BSs. Only a single macro-cell is shown for the 
sake of simplicity. 

A straightforward unifying model for heterogeneous cellular 
networks (HCNs) would consist of K spatially and spectrally 
coexisting tiers, where each tier is distinguished by its transmit 
power, BS density, and data rate, as shown in Fig. [T] For 
example, traditional BSs (tier 1) would typically have a much 
higher transmit power and lower density and offered rate than 
the lower tiers (e.g. pico and femtocells). To visualize what 
the coverage areas in such a network might look like, consider 
Figs. 2]|5 which show average power-based (equivalently aver- 
age SINR-based) coverage regions for some plausible two and 
three tier networks. Clearly, the coverage, rate, and reliability 
that mobile users experience in such networks can be expected 
to be quite different than in traditional cellular networks that 
use familiar models like the hexagonal grid. 

The objective of this paper is to provide a flexible baseline 
model for HCNs, and to show how it can be used to provide 
tractable and reasonably accurate analysis of important metrics 
like the SINR statistics, outage probability and average rate. 
Those familiar with cellular network analysis will recognize 
that this goal is fairly ambitious since such results have been 
hard to come by even for traditional cellular networks. 



A. Related Work and Motivation 

The study and design of conventional (1-tier) cellular net- 
works has often tended towards two extremes. For analysis 
and academic research, very simplistic models are typically 
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employed in order to maintain tractability, while for design 
and development (e.g. in industry) complex system-level sim- 
ulations with a very large number of parameters are generally 
used. This has made it difficult to estimate the actual gain that 
new techniques developed by researchers might provide in real 
systems. Well-known examples include multiuser detection 
(T3), multiuser MIMO (14), and BS cooperation (BJ, (16); 
all of which promised much larger gains in theory than have 
been achieved thus far in practice (T7) , (18) . 

A popular analytical model for multicell systems is the 
Wyner model JT9) , which assumes channel gains from all (usu- 
ally only 1 or 2) interfering BSs are equal and thus constant 
over the entire cell. Such a model does not distinguish between 
cell edge and interior users and in most cases does not even 
have a notion of outage since SINR is fixed and deterministic. 
It can be tuned to reasonably model average metrics in a 
system with lots of averaging, such as the CDMA uplink, but 
is not accurate in general and particularly for systems with 1 or 
2 strong interferers, like a typical OFDMA-based 4G network 
[20|. Another common approach is to consider only a small 
number of interfering cells (as few as one) and abstract the 
desired and interfering BSs to an interference channel (21) , 
[22 1 . Finally, perhaps the most popular and accepted model 
is the two-dimensional hexagonal grid model. The grid model 
is frequently used as the basis of system-level simulations but 
analysis is not generally possible |23|-|25|. However, both the 
scalability and the accuracy of grid model are questionable in 
the context of network heterogeneity (see Figs. 2]|5l 



A less accepted model is to allow the locations of the BSs 
to be drawn from a stochastic point process [26|-[28|. Such 
a model seems sensible for femtocells - which will take up 
unknown and unplanned positions - but perhaps dubious for 
the higher tiers which are centrally planned. Nevertheless, as 
Figs. [2]|5] show, the difference between randomly placed and 
actual planned locations may not be as large as expected, even 
for tier 1 macro BSs. Indeed, the recent work |29| showed 
that for a one-tier network, even with the BS locations drawn 
from a Poisson Point Process (PPP), the resulting network 
model is about as accurate as the standard grid model, when 
compared to an actual 4G network. Importantly, such a model 
allows useful mathematical tools from stochastic geometry 
to be brought to bear on the problem p0)-(32), allowing a 
tractable analytical model that is also accurate. This model 
has very recently been extended to obtain coverage results 
for femtocell networks when a typical mobile connects to its 
nearest BS |33l. 



B. Contributions and Outcomes 

The main contributions of this paper are as follows: 
General K-tier downlink model: In Section II, we define a 
tractable model for downlink multi-tier networks that captures 
many (but not all) of the most important network parameters. 
The model consists of K independent tiers of PPP distributed 
BSs, where each tier may differ in terms of the average 
transmit power, the supported data rate, and the BS density 
(the average number of BSs per unit area). The plausibility of 
the model versus planned tiers is verified through comparisons 



in Section V with an actual 4G macro-cell (1 tier) network 
with randomly placed lower tiers. 

SINR distribution, coverage probability (P c ), outage prob- 
ability (1 — P c ): Assuming (i) a mobile user connects to 
the strongest candidate BS, (ii) that the resulting Signal-to- 
Interference-plus-Noise-Ratio (SINR) is greater than 1 when 
in coverage, and (iii) Rayleigh fading, we derive an expression 
for the probability of coverage (equivalently outage) over the 
entire network under both open and closed access, which 
allows a remarkably simple closed-form in the high SINR 
regime (where interference power dominates noise power) and 
is shown to be accurate down to -4dB even under weaker 
assumptions. When all tiers have the same target SINR thresh- 
old, the coverage probability is the complementary cumulative 
distribution function (CCDF) of effective received SINR for an 
arbitrary randomly located mobile user. 
Average Data Rate: We derive the average data rate experi- 
enced by a randomly chosen mobile when it is in coverage, 
assuming interference is treated as noise but otherwise that 
the Shannon bound is achieved, i.e. the average rate in cover- 
age is E[log(l + SINR) | coverage]. This expression is readily 
computable but involves an integral so is not closed-form. 

Some interesting observations can be made from these 
results. For example, we show that when the SINR targets are 
the same for all tiers in a dense network (thermal noise power 
negligible compared to interference power), the coverage (and 
hence outage) probability does not depend upon the number of 
tiers or the density of BSs in open access, but that P c generally 
decreases with both in closed-access. This means that the trend 
towards increased density and heterogeneity and the resulting 
increase in interference need not reduce the typical SINR, 
as is commonly feared. On the contrary, aggregate network 
throughput will increase linearly with the number of BSs since 
the SINR statistics will stay the same per cell. 

We also provide the average load per tier, which is the 
average fraction of users served by the BSs belonging to a 
particular tier or equivalently the probability that a mobile 
user is served by that tier. In line with intuition, the per-tier 
load is directly proportional to the density of its BSs and their 
average transmit power, and inversely proportional to its SINR 
target. 

II. System Model 
A. Heterogeneous Cellular Network Model 

We model a HCN as a if -tier cellular network where each 
tier models the BSs of a particular class, such as those of 
femtocells or pico-cells. The BSs across tiers may differ in 
terms of the transmit power, the supported data rate and their 
spatial density. We assume that the BSs in the i-th tier are 
spatially distributed as a PPP $j of density A^, transmit at 
power Pi, and have a SINR target of fa. More precisely a 
mobile can reliably communicate with a BS x in the i-th tier 
only if its downlink SINR with respect to that BS is greater 
than fa. Thus, each tier can be uniquely defined by the tuple 
{Pii A, A^}. 

The mobiles are also modeled by an independent PPP $ m of 
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density A„J^| Without loss of generality, we conduct analysis 
on a typical mobile user located at the origin. The fading 
(power) between a BS located at point x and the typical 
mobile is denoted by h x and is assumed to be i.i.d expo- 
nential (Rayleigh fading). More complex channel distributions 
can be considered in this framework, e.g. in [29] a general 
interference fading model capable of handling any statistical 
distribution was used, and using the Fourier integral techniques 
in p4[ general fading to the selected BS can also considered. 
Both of these generalizations appear to entail significantly 
decreased tractability, and are left to future work. The standard 
path loss function is given by l(x) — \\x\\~ a , where a > 2 is 
the path loss exponent. Hence, the received power at a typical 
mobile user from a BS located at point Xi (belonging to i th 
tier) is Pih Xi ||a;i|| _a , where h Xi ~ exp(l). The resulting SINR 
expression assuming the user connects to this BS is: 

PMxiW-* 



SINR(xi) 



(1) 



where a 2 is the constant additive noise power. One of the ways 
to set the value of a 2 is according to the desired received SNR 
at the cell-edge. We will comment more on this in Section 
|Vl where we show that self-interference dominates noise in 
the typical HCNs. We assume each mobile user connects to 
its strongest BS instantaneously, i.e., the BS that offers the 
highest received SINR. Mathematically the typical node at the 
origin is in coverage if: 

maxSINR(x) > Pi, 

for some 1 < i < K. An assumption that greatly simplifies 
the analysis is that ft > 1 ( dB^] The following Lemma 
shows that under this assumption, at most one BS in the entire 
network can provide SINR greater than the required threshold. 
Although some users in commercial cellular networks indeed 
have operating SINR below dB, they are in a distinct minority 
(cell edge users) and in Section V we show numerically 
that this model holds very accurately at least to —4 dB, 
which covers cell edge users as well. The following Lemma 
characterizes the number of potential BSs that a mobile can 
connect to and will be used in the later sections. 

Lemma 1. Given positive real numbers {oi, ai . . . a n }, which 
correspond to the received power from each BS at the typical 
mobile user and defining Ci = ^ a ^ +g . 2 , which corresponds 

to the SINR of the i BS, at most m Ci's can be greater than 
1 jm for any positive integer m. 

Proof: See Appendix. ■ 

B. Coverage Regions 

Before going into the analysis and main results, it may 
be helpful to first build some intuition about the proposed 
model, and its resulting coverage regions. The illustrative 

'Since the SINR distribution is derived as a spatial average, the distribution 
of users is implicitly assumed to be homogeneous. 

2 This assumption is relaxed in [35] where we compute the coverage 
probability for general ft. This also enabled us to compute the ergodic rate 
for a typical mobile user. 




Fig. 2. Coverage regions in a two-tier network as per the model used in 
this paper. Both macro (large circles) and femto (small dark squares) BSs are 
distributed as independent PPPs with Pi = IOOOP2 and A2 = 5Ai. 




Fig. 3. Coverage regions in a two-tier network where Macro (tier- 1 ) BS 
locations (large circles) correspond to actual 4G deployment. Femto BSs 
(small dark squares) are distributed as a PPP (Pi = IOOOP2 and A2 = 5Ai). 



HCN coverage regions can be visually plotted in two steps, 

M 



resulting in Figs. 



First, we randomly place K different 
types of BSs on a 2-D plane according to the aforementioned 
independent PPPs. Ignoring fading, the space is then fully 
tessellated following the maximum SINR connectivity model, 
which is equivalent to maximum SIR and maximum power 
connectivity models in the absence of fading. Please note that 
in reality the cell boundaries are not as well defined as shown 
in these coverage regions due to fading. Therefore, these plots 
can be perceived as the average coverage footprints over a 
period of time so that the effect of fading is averaged out. 
Due to the differences in the transmit powers over the tiers, 
these average coverage plots do not correspond to a standard 
Voronoi tessellation (also called a Dirichlet tessellation) |36|. 
Instead, they closely resemble a circular Dirichlet tessellation, 
also called a multiplicatively weighted Voronoi diagram (37). 
The coverage regions for a two-tier network - for example 
comprising macro and femtocells - are depicted in Figs. [2] 
and [3] for two cases: 1) the macro-cell BSs are distributed 
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Fig. 4. Close-up view of coverage regions in a three-tier network. All the 
tiers, i.e., tier-1 macro (large circles), tier-2 pico (light triangles) and tier-3 
femto (small dark squares), are modeled as independent PPPs. Pi = IOOP2 = 
IOOOP3, A3 = 4A 2 = 8Aj. 




Fig. 5. Coverage regions in a three-tier network where macro BS locations 
(large circles) now correspond to actual 4G deployment. Other parameters are 
same as Fig. [4] 



according to PPP (our model), and 2) the macro-cell BSs 
correspond to an actual 4G deployment over a relatively flat 
urban region. The femtocells are distributed according to an 
independent PPP in both cases. Qualitatively, the coverage 
regions are quite similar in the two cases. 

In Figs. [4] and |5J the coverage regions are now shown 
with an additional pico-cell tier. As is the case in the actual 
networks, we assume that the macro-cells have the highest 
and the femtocells have the lowest transmit power, with 
pico-cells somewhere in between. For example, in LTE |38|, 
typical values are on the order of 50W, .2W, and 2W, re- 
spectively. Therefore, femtocell coverage regions are usually 
much smaller than the other two tiers, particularly when they 
are nearby a higher power BS. Similarly, we observe that 
the coverage footprint of pico-cells increases when they are 
farther from the macro BSs. These observations highlight the 
particularly important role of smaller cells where macrocell 
coverage is poor. 



C. Applicability of the Model 

The model is applicable both to non-orthogonal 
and orthogonal (TDM A, OFDM A) cellular networks. The 
analysis is for a single frequency band and assumes that all 
BSs are transmitting continuously in all time slots at constant 
power, although if some fraction / of time slots were not 
used (at random), then the resulting density of interfering BSs 
would simply be (1 — /)A and the analysis could be extended. 
In OFDMA-based networks, it is desirable to move strongly 
interfering neighbors or tiers to orthogonal resources in time 
and/or frequency and so the coverage can be improved. Sim- 
ilarly, additional enhancements like opportunistic scheduling 
or multiple antenna communication should increase coverage 
and/or rate and this framework could be extended to indicate 
the gains of different approaches. Although we do not explic- 
itly consider antenna sectoring, it can be easily incorporated in 
the current model if sectoring is done randomly. If the beam 
is partitioned into n equal sectors, the density of interfering 
BSs reduces by a factor of n because the probability that the 
beam of any BS would point towards a randomly chosen BS 
is 1/n. Cellular engineers will note that further details are 
missing from this model. In addition to shadowing, we do not 
consider frequency reuse, power control, or any other form of 
interference management, leaving these to future extensions. 
In short, this is a baseline tractable model for HCNs. 

III. Coverage Probability and Average Load per 

Tier 

A typical mobile user is said to be in coverage if it is able to 
connect to at least one BS with SINR above its threshold. In the 
case when all the tiers have same SINR threshold /3 > 1, cov- 
erage probability is precisely the complementary cumulative 
distribution function (CCDF) of the effective received SINR, 
outage being the CDF, i.e., 1— CCDF. With this understanding, 
we now derive the probability of coverage for a randomly 
located mobile user both for open and closed access networks 
(defined below). Using these results, we also derive a measure 
of average load per tier in terms of the fraction of users served 
by each tier. 

A. Open Access 

We first assume the open access strategy where a typical 
mobile user is allowed to connect to any tier without any 
restriction. Under the current system model, this strategy 
reduces to choosing the strongest BS, i.e., the one that delivers 
the maximum received SINR. 

1 ) Coverage Probability: The main result for the probabil- 
ity of coverage in open access networks is given by Theorem[T| 

Theorem 1 (General case). When fa > I, the coverage 
probability for a typical randomly located mobile user in open 

3 For CDMA networks, although the received SINR is generally much 
smaller than 1, the post-despreading SINR, which is what determines cov- 
erage/outage, is often greater than or at least close to 1, so our model and 
Lemma 1 are still reasonable if the interference term is divided by a spreading 
factor M. 
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access is 



K . 

Pc({A4, {Pi}, {Pi}) = E A * / ex P ( - 



2/c 



if 



II^H 2 J] A m P^ Q exp 



ft<r 2 



||a; l |r)da; l , 



(2) 



w/zere C(a) = 2tt 2 csc(^)a~ 1 . 

Proof: See Appendix. ■ 
Theorem [T] gives a simple and fairly general expression for 
coverage probability. For better understanding of the proof, we 
now provide a brief description of the main steps. First recall 
that a mobile user is in coverage if it is able to connect to at 
least one BS with SINR above its threshold. Now assuming, 
ft > 1, V i, we know from Lemma[T]that a mobile can connect 
to at most one BS. Therefore, P c can now be defined as the sum 
of the probabilities that each BS connects to the mobile (with 
the understanding that all the events are mutually exclusive and 
at most one of them happens at any time). This leads to a sum 
of probabilities over PPP, which can be converted to a simple 
integral of Laplace transform of cumulative interference using 
Campbell-Mecke Theorem pO) . A closed form expression for 
the Laplace transform can be evaluated in two main steps. 
Firstly, the nature of interference function (sum over PPP) 
leads to a product form for its Laplace transform. Using 
probability generating functional (PGFL) of PPP |30) and the 
fact that fading power is exponentially distributed, we arrive 
at the closed form expression for Laplace transform which 
directly leads to the final result of the Theorem. This result can 
be simplified further for the interference-limited case, where it 
reduces to a remarkably simple closed-form expression given 
by Corollary [T] 

Corollary 1 (No-noise). In an interference limited network, 
i.e., when self-interference dominates thermal noise, the cov- 
erage probability of a typical mobile user simplifies to 



P c ({A 4 },{ft},{PJ) = 



2/a 



ft > 1. 



Proof: Follows from Theorem [I] with a 2 = 0. ■ 
The simplicity of this result leads to some important obser- 
vations. Firstly, setting K = 1 leads to the single-tier case, 
where the coverage probability is given by: 



Pc(A,/?,P) = 



(3) 



C(a)f3 2 / a ' 

From Q, we note that the P c in an interference-limited single- 
tier network is independent of the density of the BSs, and 
is solely dependent upon the target Signal-to-Interference- 
power-Ratio (SIrQ This is consistent with |29|, where a 
similar observation was made for a single-tier network using 
nearest neighbor connectivity model. The intuition behind this 
observation is that the change in the density of BSs leads to 

4 When the system is interference-limited, SINR and SIR can be used 
interchangeably since thermal noise is negligible compared to the interference 
power. However, for concreteness we will henceforth use SIR when assuming 
interference-limited network. 



the change in the received and interference powers with the 
same factor and hence the effects cancel. 

From Corollary [T] it also follows that, if ft = ft V i, 
in an interference-limited network then P c ({Ai}, ft {Pi}) = 
c(a)p 2 / a • This * s P erna P s an unexpected result since it states 
that the coverage probability is not affected by the number 
of tiers or their relative densities and transmit powers in an 
interference-limited network. In fact, it is exactly the same as 
that of the single-tier case. Therefore, more BSs can be added 
in any tier without affecting the coverage and hence the net 
network capacity can be increased linearly with the number of 
BSs. The intuition behind this result is that the decision of a 
mobile user to connect to a BS depends solely on the received 
SIR from that BS and a common target SIR, unlike the general 
case where it also depends upon the tier to which the BS 
belongs. Thus, the mobile user does not differentiate between 
the tiers when the SIR thresholds are the same. Surprisingly, 
this leads to a situation similar to the one discussed for the 
single-tier case above, where the change in the received power 
due to the change in the density or the transmit power of BSs 
of some tier is equalized by the change in the interference 
power. Significantly, this implies that the interference from 
smaller cells, such as femtos and picos, need not decrease 
network performance in open access networks. 

2) Average Load per Tier: The average load on each tier is 
defined as the average fraction of users in coverage served by 
that tier. This can also be interpreted as the average fraction of 
time for which each mobile is connected to the BSs belonging 
to a particular tier. The main result for the average load per 
tier in open access is given by Proposition [T] 

Proposition 1. The average fraction of users served by j th 
tier (also the average load on j th tier) in open access is 



A, 



Pc({AJ,{ft},{P,}) 

K 

INI 2 ^ A m P 2 Mexp 



exp 



(-C(a) 



2/c 



fta 2 
Pi 



\xA\ a )dx 



Proof: See Appendix. ■ 
In an interference-limited scenario, this result reduces to a 
simple closed form expression, which is given by the following 
Corollary. 

Corollary 2. When noise is neglected, i.e., a 2 = 0, 



*3 = 



j 



\P- ,a pf /a ' 



From Corollary [2j we observe that the load on each tier 



2/a o-2/a 



In line 



is directly proportional to the quantity XjPj ftj 
with intuition, a tier will serve more users if it has a higher 
BS density or higher transmit power or a lower SIR threshold. 
When the thresholds of all tiers are equal to j3 and the transmit 
powers of all BSs equal to P, the average load on each tier is 

— \ 

. Hence, as expected the average load on each 



— 22 



tier is directly proportional to the density of its BSs. 
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B. Closed Access 

Under closed access, also known as a closed subscriber 
group, a mobile user is allowed to connect to only a subset 
of tiers and the rest of the tiers act purely as interferers. The 
motivation for closed access particularly applies to privately 
owned infrastructure, such as femtocells or perhaps custom 
picocells mounted on a company's roof to improve service 
to their staff. The desirable aspects of closed access can 
include protection of finite backhaul capacity, security, and the 
reduction in the frequency of handoffs experienced by mobile 
users and the associated overhead required. In the context of 
our model, closed access means that if the strongest BS lies 
in the restricted tier, it by definition leads to an outage event 
irrespective of the received SINR associated with that BS. Fur- 
thermore, since closed access is a constraint on connectivity, 
it would naturally lead to reduced coverage probability. This 
intuition is verified in the following discussion. 

1) Coverage Probability: The main result of coverage 
probability in closed access networks is given by Lemma [2] 

Lemma 2. When a typical mobile user is allowed to connect 
to only a subset B C {1,2,,..., K}, the coverage probability 
for closed access is 



Pc(M, {A}, {Pi}) =J2 Xi [ ex P ( - c ^ (§) 



2j, 



K 



|N| 2 £A m P£/ Q )exp(- 



m— 1 



Pi 



01 Ida*. 



(4) 



Proof: The coverage probability is 



Pc({Ai},{A},{Pi})=P |J SIHR(xi)>A 

-E E E [l(SINR(^)>A)], 

where (a) again follows from Lemma [T] under the assumption 
that Pi > 1. Following the same steps as the proof of Theorem 
[T] we arrive at the final result. ■ 
The following corollary specializes from Lemma [2] to 
interference-limited HCNs. 



Corollary 3. When a 2 = 0, 



PcttAiMAM^}) 



C(a) 



El^iPi 



2/a 



If the threshold of each tier to be same (and equal to (3) 
and the transmit power of each tier to be same (and equal 



£ 4e 



'. So, if the 



to P), the coverage probability is c ^^ 2/a ^ K A 
thresholds and transmit powers of all the tiers are same, closed 
access has a lower coverage than open access by a factor of 

Ef =1 A, ' 

2) Average Load per Tier: The main result for the average 
load per tier under closed access is given by Proposition 
The proof directly follows from the proof of Proposition 
with the understanding that the coverage event would now be 
defined by only the "allowed" tiers. 



Proposition 2. The average fraction of users in coverage 
served by j th tier (also the average load on j tier) in closed 
access is 



N, 



Pc({A 4 },{ft},{P,}) 




otherwise. 



(5) 



where P c ({Aj}, {A}, {Pj}) is the coverage probability under 
closed access given by Lemma [2] and 

5 = [ e-^r^lNI 2 ^!^" 7 " 



e j ax. 



K 2 



The corresponding result for the interference-limited networks 
is 







ieB, 

otherwise. 



(6) 



IV. Average Rate 

In this section, we derive the average rate R achievable 
by a random mobile user when it is in coverage both for 
the open and closed access strategies. It is worth noting that 
since the rate is computed conditioned on the mobile being 
in coverage, it is not the same as the classic ergodic rate 
E[R}. The motivation behind considering this metric is that 
given the coverage/outage information, the service providers 
are interested in knowing the average rate they can provide to 
the users that are in coverage. 

A. Open Access 

The main result for the average rate in open access is given 
in Theorem [2] In this section, for notational simplicity, we 
restrict our attention to the case of a 2 = 0. However, the 
results can be extended to the general case with noise in a 
straightforward manner. 

Theorem 2. The average rate achievable by a randomly 
chosen mobile in open access when it is in coverage is 

-ft — log (1 + Prnin) H _ ^ . _ 9/ oln , (7) 



where 



max(Pi,x) 2 / a 



Ax, 



and f3 min = min{/3i, f3 2 , Pk}- 



Proof: See Appendix. ■ 
Thus we observe that the average rate expression involves 
only a single integral which can be easily evaluated numeri- 
cally. 

Corollary 4. Using the same threshold ft for all tiers, the 
average rate achievable by a randomly chosen mobile that is 
in coverage in open access is: 

R = log(l + f3)+f3 2 ^ a A(a,f3,f3). (8) 

The above result shows that the average rate is independent 
of the density of BSs of each tier when the SIR thresholds are 
same for all the tiers. This is expected because the distribution 
of max SIR does not depend upon the density of BSs in this 
case (follows from Theorem [TJ. 
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B. Closed Access 

The average rate R c achievable by a randomly chosen 
mobile under closed access (assuming it is under coverage) 
can be expressed as: 



E 



log 1 



max (SIR(sg)) 



|J |J (SIR(x) > ft 



(9) 

Following the same steps as in proof of Theorem [2] we arrive 
at the following Proposition. 

Proposition 3. Assuming a mobile user is allowed to connect 
to only a subset B of the K tiers, the average rate ( assuming 
mobile is under coverage) can be expressed as: 



R c = log (1 + ft, 



heb 



XiPf /a p- 2/a 



(10) 



where ft 



min{ft}. 

lEB 



Corollary 5. Assuming the threshold of each tier is the same 
and equal to ft tlie average rate achievable by a randomly 
chosen mobile in coverage under closed access is 

R c = log(l + ft + 2 ' a A{a, ft ft. (11) 

From Corollaries |4] and [5] we observe that the average rate 
(R) of the mobile while it is in coverage is not affected by 
access control when the thresholds are the same for all tiers. 
However, since the coverage probability is lower in case of 
closed access, it would naturally lead to a lower ergodic rate 
as compared to the open access networks. Interested readers 
can refer to (35) for the derivation of ergodic rate in this 
framework. 

V. Numerical Results 



Simulation: PPP 
Simulation: Actual BSs 
Simulation: Grid 
Analysis 




2 4 

SIR Threshold (p ) 



Fig. 7. Coverage probability in a two-tier HCN (K 
100P 2 , A 2 = 2Ai, fa = 1 dB, No noise). 



2, a = 3, Pi 



Simulation: PPP 
Simulation: Actual BSs 
Simulation: Grid 
Analysis 




2 4 

SIR Threshold (p ) 




-6 -4 -2 2 4 6 8 10 

SIR Threshold (p ) 



Fig. 6. Coverage probability in a two-tier HCN with and without thermal 
noise (K = 2, Pi = 25P 2 , A 2 = 5Ai, fa = 1 dB, SNR edge = dB). 

Most of the analytical results presented in this paper are 
fairly self-explanatory and do not require a separate numerical 
interpretation. Therefore, to avoid repetition, we will present 
only non-obvious trends and validation of the model in this 
section. 



Fig. 8. Coverage probability in a two-tier HCN (K = 2, a = 3.2, Pi = 
1000P 2 , A 2 = 4Ai, fa = 1 dB, No noise). 

A. Effect of Thermal Noise 

We first study the effect of thermal noise on the coverage 
probability by considering a typical two-tier network consist- 
ing of macro-cells overlaid with pico-cells. To set the noise 
power, we use the following notion of cell-edge users in this 
example. Defining the distance of the the nearest macro BS 
to the typical mobile user to be d and the underlying random 
variable to be D, the mobile user is said to be on the cell 
edge if P(D < d) > P e dge, where P e d ge is set to 0.9 for 
this illustration. For PPP(A), P(D < d) = 1 - cxp(A7rd 2 ), 

giving d > V ' n(1 ~A Pe,ia ^ For a desired ed 8 

e-user SNR, say 

SNR e rf S e, a 2 can be approximated as a 1 w ^ R "^ ge , where 
d e dge is the limiting value of d evaluated above. Under this 
setup, we present the coverage probability for various values 
of a in Fig. [6] By comparing these results with the no-noise 
case, we note that the typical HCNs are interference limited 
and hence thermal noise has a very limited effect on coverage 
probability. Therefore, we will ignore noise in the rest of this 
section. 
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B. Validity of PPP Model and /3 > 1 Assumption 

While a random PPP model is probably the best that can be 
hoped for in modeling "unplanned" tiers, such as femtocells, 
its accuracy in modeling "planned" BS locations, such as those 
of macro-cells, is open to question. Therefore we verify the 
PPP assumption for macro-cells from a coverage probability 
perspective by considering a two-tier network in three different 
scenarios: 1) the macro-cell BSs are distributed according to 
PPP (our model), 2) the macro-cell BSs correspond to an 
actual 4G deployment, and 3) macro-cell BSs are distributed 
according to hexagonal grid model. The second tier is modeled 
as an independent PPP in all three cases. As shown in 
Figs. [7J and [8] the actual coverage probability lies between 
the coverage probabilities of the PPP and grid model. This 
is because the likelihood of a dominant interferer is highest 
for the PPP and lowest for the grid model. This comparison 
shows that the PPP assumption is nearly as accurate as the 
grid model in the case of macro-cells, with the PPP providing 
a lower bound and grid model providing an upper bound to 
the actual coverage probability. 

We now focus on the /3 > 1 assumption by comparing 
the theoretical and simulated results for coverage probability 
in Figs. [JJ and [8] As expected, the simulated and analytical 
results match reasonably well for ft > 1 but interestingly, the 
theoretical results also provide a tight upper bound to the exact 
solution even until about ft = —4 dB (sa .4). Therefore, the 
analytical results also cover typical cell edge users. The same 
trend is observed in the case of average rate results presented 
in Fig. [9] which are also accurate down to about —4 dB target- 
SIR. 




SIR Threshold (P) 



Fig. 9. Average rate while mobile is in coverage (K = 2, a = 3, Pi 
1000P2, -^2 = 2Ai, fii = /?2 = no noise, open access). 



it would be useful to include scheduling, resource allocation 
and/or various forms of frequency reuse (such as fractional 
frequency reuse) in this setup. At the network level, further 
BS cooperation techniques could use the framework. 

The random spatial model used in this paper could likely 
be further improved by incorporating a point process that 
models repulsion or minimum separation distance between 
BSs, such as determinantal and Matern processes (30), [3~T) , 
respectively. Often picocells or enterprise (operator deployed) 
femtocells are clustered in high-demand areas, so a Poisson 
cluster process [39| might be useful to model such scenarios. 
It would also be very helpful to consider further actual 
deployments from the one considered in this paper (fiat, large, 
low density urban), and see which of these general point 
processes best model different deployment scenarios such as 
dense urban (New York, Tokyo) vs. sprawling (Los Angeles, 
Sydney), fiat relatively uniform cities (Paris, London) vs. 
mountainous/coastal cities (Rio de Janiero, Hong Kong), and 
so on. Finally, the user distribution was implicitly assumed 
in this paper to be homogeneous, with all BSs actively 
transmitting at all times. An important extension is to con- 
sider non-homogeneous user distributions and realistic traffic 
models to better understand optimal cell association policies 
in heterogeneous cellular networks. An initial characterization 
of biased cell association policies can be found in pO) . 

Appendix 

A. Proof of Lemma [7] 

Since SINR < SIR for all the BSs, defining b % = — 

l^j^i a 3 

as the SIR corresponding to the i BS, it suffices to show 
that at most m fej's can be greater than 1/m for any positive 
integer m. This is shown below. 



h = 



E 

bi 
l + b 

n 

E 



E 



j 3 



= 1. 



(12) 



We first prove the result for m = 1 (by contradiction) and then 
show that it can be trivially extended to the case of general in. 
We first observe that ( fT2] > is satisfied if only one of the frj's is 
greater than 1. Now assume that two 6;'s are greater than one 
and without loss of generality, assume that they are bi and 6 2 . 
This implies 1/bx and l/6 2 € (0,1). Therefore, and 
57^6(1/2,1). Thus, 



VI. Conclusions 

We have provided a new tractable model for if-tier down- 
link HCNs and an associated analysis procedure that gives 
simple mathematical expressions for the most important per- 
formance metrics. The possible extensions of this work are 
numerous, and could include physical layer technologies like 
multiple antennas, spread spectrum, power control, interfer- 
ence cancelation or interference alignment. For the MAC layer, 



1/& 



1/6, + 1 f-^ 1/h + 1 

l — l 

n 1 

>1 + §I7^TT 



(13) 



which is in contradiction with ( fT2| ). Since (V2\ does not even 
hold for two 6, 's greater than one, it proves that the only one 
of the 6j's can be greater than one. Similarly for the case of 
general m, it is easy to observe that ( fT2] i is trivially satisfied if 
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at most m of the bi's are greater than 1/m. Now assume that 
m + 1 bi's are greater than 1/m and without loss of generality, 
assume that they are b\, b 2 , . . . , 6 m +i- Proceeding as in ( pj) , 



n 1 
^ + 1 

2—1 



> 1 



i=m+2 ' 



(14) 



which is in contradiction to | [T2) . Therefore, at most m 6j's 
can be greater than 1/m. ■ 

B. Proof of Theorem [JJ 

For notational simplicity, denote the set {1, 2, . . . K} by /C. 
The coverage probability in a if -tier network under maximum 
SINR connectivity model can be derived as follows: 

Pc({AJ,{ft},{P}) 

|J SIMftfo) > ft 



E 



1 (J SINRfcj) > ft 



(a) 



A' 



= E E E [1 (SINR(^) > ft)] 
i=l 

ffl^/ pf P ^^)>ft) dXi 

^ JUL* \ Ixi + V 1 

K e / a 

i= i 



p/(x 4 ) 



(15) 



where (a) follows from Lemma [T] under the assumption that 
ft > 1 V i, (b) follows from Campbell Mecke Theorem |30|, 
and (c) follows from the fact that the channel gains are 
assumed to be Rayleigh distributed. Here Cj x , (.) is the 
Laplace transform of the cumulative interference from all the 
tiers when the randomly chosen mobile user is being served 
by the i th tier. Since the point processes are stationary, the 
interference does not depend on the location Xi. Therefore, 
we denote Ci x by Ci i , which is given by 



K 



C u (s) = ]Je Ii 

3=1 



Y\ exp (sPjh Xj l{xj)) 



Using the independence of the fading random variables Cj i (s) 
equals 



K 

3=1 



K 



3=1 
K 



n 

11 1 + sP 



xj£<bj/xi 



+ sPjl{Xj) 



= TT exp ( -Xi / 



X + sPjWxjW- 01 



dx. 



( = } n cx P (-27rA i ( s P J -) 2/Q / r 
■ =1 V Jo Jo 



e (-*(l+^)) d t drj 
(16) 



where (a) follows from the Rayleigh fading assumption (i.e., 
h ~ exp(l)), (6) follows from probability generating func- 
tional (PGFL) of PPP 1 30 1 and, (c) results from algebraic ma- 
nipulation after converting from Cartesian to polar coordinates 
Using some properties of Gamma function, ( p"6| ) can be further 
simplified to 



K 



a/a i 



i=l 



where C(a) = 



2ir 2 csc(^) 



probability P C ({AJ, {ft}, {P}) is 
K 



Using ( fl"5j ) and ( fT7| ) the coverage 



E* 



e -^J c( Q )|| :Cs |I^Ef„ =1 A m P^° e ^|K|l d;E . 



which completes the proof. ■ 

C. Proof of Proposition [JJ 

Let P n 6 R 2 denote an increasing sequence of convex sets 
with B n C B n+1 and linin^oo |P„| = oo. For this proof, 
we denote STR Xrn (%b) as the received SIR when the mobile is 
located at x m ^ connects to BS located at Xf,- Please recall 
that the subscript is dropped and SIR is denoted as SIR(xfe) 
when the mobile user is located at the origin. The average 
fraction of users served by the j th tier can now be expressed 
as: 



Nj = hm -, r > 



1 (J SINR Xm (>,,•) > ft I (J (SINR Xm (a: i ) > ft) 



(a) 



(&) 



U SINRfo) 



> 



(J (SINR^) > ft) 



U SINR(^) > ft, U (SINR(^) > ft) 



^U eM *,( SINR (^) >ft) 
:Ux j6 », SINR(^) > ft 



u. 



£K,Xi£<&i 



(sim(xi) > ft 



(18) 



where (a) follows from the stationarity and the ergodicity 
of PPP 1301. P !o denotes the reduced Palm distribution of a 



PPP and (b) follows from the Slivinak's theorem [30], |41| 
and Bayes rule. Noting thatP \[j Xj ^. SINR^j) > ftj is me 
probability of coverage with a single tier j, the result follows 
from Theorem Q] ■ 
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D. Proof of Theorem [2] 

Denoting the coverage event U»=i Uxe^i (SIR(x) > Pi) by 
C({/3j}), the average rate achievable by a randomly chosen 
mobile user when it is under coverage can be expressed as: 



R = E 



log 1+ max (SIR(x)) C({/3J) 



(19) 



We first derive the conditional complementary cumulative 
density function (CCDF) of max^y $i (SIR(x)) as follows: 

max (SIR(x)) > T C({/3J) 
(a) P (max^y g, (SIR(») > T, C({ft») 

P(C({A») 

(6) P(C({T»,C({A») 
P(C({ft})) ' 
_ P(C({max(T,ft)})) 



P(C({A})) ' 



otherwise 



(20) 



where (a) follows from Bayes' theorem, (6) follows from 
Lemma [T| under the assumption ft > 1 V i, (c) follows from 
Theorem 111 and /3 TO j n denotes rnin{/?i, /?2, ■ ■ ■ , /3k}. 

Denoting random variable max^gy (SIR(x)) by X, R can 
be evaluated as follows: 



R 



log(l + x)f x (x\C({P i }))dx, 



OO FX 



-f x {x | C({ft})) dy dx, 



(«) 



a;=0 «/j/=0 1 + V 

/ / / x (x | C({&})) dx) — — dy, 
V(X>y\C({(3 % })) 



l + y 



-dy, 



(21) 



where (a) follows from changing the order of integration. Now 
we substitute ( |2"0] > in ( f2"T| to get the average rate as: 



R 



l + y 



dy 



1 



K 



= l0g(l + /3 m m) + 
1 



a/Q /°° max(ft,i)- 2 / a 



1 + x 



This completes the proof. 



j2/«fl-2/ 



maxfft.j:)- 2 / 



1 + x 



dx 



(22) 
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